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A BRIEF ANALOG INVESTIGATION OF INERTIA COUPLING IN 
ROLLING MANEUVERS OF AN AIRPLANE CONFIGURATION 
USING A VARIABLE -INCIDENCE WING AS 
THE LONGITUDINAL CONTROL 
By Clarence L. Gillis 


SUMMARY 


An analog computer study was made in order to investigate the iner- 
tia coupling in rolling maneuvers of an airplane configuration utilizing 
a variable-incidence wing for longitudinal control. A five -degree -of - 
freedom system of eqtiations was used, and only one set of flight condi- 
tions was included in this study. The resvilts, when compared to those 
for a conventional tail -control airplane, indicated significant advan- 
tages for the variable -incidence -wing type of control in reducing unde- 
sired dynamic effects during rolling maneuvers . 


INTRODUCTION 


The problem of inertia coupling of the longitudinal and lateral 
motions of airplanes in rolling maneuvers has been extensively investi- 
gated, both experimentally and analytically, and various methods have 
been proposed for eliminating roll coupling or at least reducing the 
magnitude of the coupling to acceptable values. (See refs. 1 to 7*) 
Reference 7 contains a good discussion of the effects of the individual 
coupling terms and of the interaction of these coupling terms to produce 
the complex dynamic response during rolling maneuvers. The investigation 
described herein concerns a proposed method of operation of airplanes 
which might lead to reduced inertia coupling as well as to some other 
possible advantages. The method of operation consists of obtaining nor- 
mal acceleration by varying the wing incidence while keeping the fuselage 
at a relatively constant angle of attack, as is done for some missiles. 
Other studies (refs. 5 a-nd 7 ) have shown that one of the most important 
effects causing coupled longitudinal and lateral motions is a term in 
the side-force equation proportional to the product of angle of attack 
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and rolling velocity. By permitting the fuselage^ and thus the princi- 
pal axis, to remain at a small and relatively constant angle of attack, 
the variable-incidence-wing type of longitudinal control might thus 
alleviate the inertia coupling problem. 

The present investigation is preliminary because the investigation 
considers only one configuration, which may not be optimum as far as 
variable-incidence operation is concerned, and because the calculations 
are made for only one set of flight conditions aind for one type of 
rolling -pullout maneuver. Calculations to simulate airplane motions 
were made on an electronic analog computer at the Langley Aeronautical 
Laboratory, using a five -degree -of -freedom system of equations, with 
forward velocity considered constant. The flight conditions and mass 
characteristics used for the calculations represented a free-flight 
rocket-propelled model airplane. The application of results to a full- 
scale airplane is discussed. 


SYMBOLS 


*^m 


lift coefficient. 


Lift 

q-S 


pitching -moment coefficient. 


Pitching moment 
q'Sc 


rolling -moment coefficient, moment 

q'Sb 


Cn 


C 


Y 


yawing -moment coefficient, 
lateral-force coefficient. 


Yawing moment 
q'Sb 

Lateral force 
q'S 


b 

c 




''xz 


wing span, ft 

wing mean aerodynamic chord, ft 

moments of inertia about X-, Y-, and Z-axes (fig. l), 

slug-ft2 

product of inertia 


k^,ky,k2 


radii of gyration about X-, Y-, and Z-axes, ft 
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M 

m 

P 

Pcc 

q 

r 

5 
t 
V 
W 
a 

P 

6 

P 

P 


Mach number 
mass^ slugs 

rolling angular velocity, radians/sec 
atmospheric pressure, Ib/sq. ft 

pitching angular velocity, radians/sec 

dynamic pressure, Ib/sq ft 

yawing velocity, radians/sec 

wing area, sq ft 

time, sec 

velocity, ft /sec 

weight, lb 

angle of attack, radians 
angle of sideslip, radians 

longitudinal control deflection (positive with trailing 
edge down) , radians 

differential deflection of two ailerons (positive in a 
direction to give positive roll) , radians 

atmospheric density, slug/cu ft 

relative density factor, m/pSc 


Subscripts: 

max 

s 

ss 

w 

c 


maximum 
stabilizer 
steady state 
wing 

conventional configuration 


v 


var iab le - inc idenc e -wing c onf igurat ion 
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The quantities a, d, and 6a used as subscripts to coeffi- 
cients indicate the partial derivative of the coefficient with respect 
to the quantity, such as 



The quantities a, 3, p, q, and r used as subscripts to coefficients 
indicate the partial derivative of the coefficient with respect to the 
nondimens ional quantities dc/2V, 3b/2V, pb/2V, qc/2V, and rb/2V 
such as 



and Ct 
'•P 




EQUATIONS OF MOTION 


The equations of motion were written for the body axes system illus- 
trated in figure 1. These axes are considered fixed in the fuselage for 
all calculations presented herein. The eq\mtions as set up for the ana- 
log computer were 


d = q - pp - 


CL'S 

mV 






( 1 ) 


3 = 


-r + ap + 


q'S 

mV 




( 2 ) 


P = 


% 


qr 


^XZ N q', 

+ — r + pq + — 
ix V ^ I 


q'Sb 


("l6> * "ipP) 


q'Sb^ 

2VIx 



( 5 ) 




q = 


P, , ^ . pa) , (S" - ^ fC: 




cc + Cjngb 


(M 
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r 



+ 


CnpP + 



a'SV 

2VIz 


2 



( 5 ) 


_r + Cn-P + Cn^p 


All gravity terms have been omitted from the above eqtiations as well as 
certain aerodynamic terms such as Cj„(p) and Cy^. For the flight con- 
ditions used herein the gravity terms are small and approximately equal 
for the two configurations, and the aerodynamic terms omitted are assumed 
to be negligible. In the calculations it was ass\jmed that p = -r in 
equation (5), thus a slight si2ig)lification of the last term is possible. 


The calculations were made for the config\nration shown in figure 2 . 
This configuration was chosen from considerations of expediency, and no 
studies were made to derive an optimum configuration when a variable- 
incidence wing was enployed. Aerodynamic data on this and a similar 
configuration were available from references 8 and 9» 

The calculations were made for a Mach nximber of 2 . 0 . The flight 
conditions used were those existing on the rocket -powered model of refer- 
ence 8 at this Mach number. These conditions, at low altitudes (see 
table l), are dynamically similar to a representative f\lLl-scale airplane 
at altitudes in the range of 50,000 to 50,000 feet, as is discussed in 
the section on "Application to Fuill-Scale Airplane." 

Calculations were made for both a conventional configuration, that 
is, a configvoration with a fixed wing and an all -movable stabilizer, and 
for a configuration with an all -movable wing and a fixed stabilizer. The 
aerodynamic characteristics used for the two configurations are given in 
table II. In estimating and 0 ^^ for the variable -incidence -wing 

airplane the loss in lift across the fuselage and the effect of the down- 
wash on the tail from the deflected wing were accounted for. No changes 
in the derivatives with angle of attack were considered. The calcuila- 
tions included two values of as indicated in table II. 

In order to keep the fuselage at a relatively fixed angle of attack 
during maneuvers, the term ^mg for the variable -incidence wing control 

ideally must be zero. For any given set of flight conditions can 


CALCULATIONS 


6 


MCA RM L57F18 


usually "be kept at zero "by making fairly small changes in configuration; 
for example, in the present case can he made zero while keeping 

the same as for the conventional configuration by moving the wing 

forward 7 percent of the mean aerodynamic chord and approximately 
doubling the horizontal -tail area, as indicated by the dotted lines in 
figure 2. Further comments on this subject are included subsequently. 

The maneuver chosen for the calculations was a pull-up followed by 
a roll. The airplane was assumed to be in level flight initially, and 
an abrupt stabilizer or wing deflection was applied to produce a change 
in steady-state lift coefficient of 0.2. Aileron deflections were then 
applied at two different times dinging the maneuver. In one set of runs 
the transient motion resulting from the pull-up was allowed to damp to 
a small magnitude, and the aileron deflection was then applied abruptly. 
In the other set of runs the aileron deflection was applied abruptly at 
or near the first peak in normal acceleration (or angle of attack) fol- 
lowing the control motion. In both cases the ailerons remained deflected 
for a time interval sufficient to give approximately a 150^ to l80^ roll. 
The ailerons were then moved abruptly back to neutral, and following 
another time interval to permit damping of the resulting transient the 
stabilizer was also returned to neutral. Some studies have indicated 
that larger transient motions may be encountered in recovering from a 
maneuver than during the maneuver. Several different aileron deflec- 
tions were used to produce rolling velocities from low values up to 
values near those corresponding to resonance with the yawing motion 
(20 radians/sec). The natural frequency in pitch was 25 radians /second. 
The rolling divergence at resonance was not investigated in this study. 
Almost all of the operational rolling motions of an airplane occur at 
roll rates less than the resonance rate, and even under these conditions 
severe excursions in pitch and yaw can occur and lead to violent motions 
and excessive aircraft loads (ref. k ) . Calculations for the variable- 
incidence -wing airplane included several other conditions, but these 
conditions are more conveniently discussed with the results where their 
p-urpose becomes clear. 


RESULTS 


A total of about 60 analog computer rims were made. For illustra- 
tive purposes 25 of these rims are shown in figures ^ to ^ as time his- 
tories of lateral and longitudinal motions . Summary plots of some of 
the significant results are presented in figures 8 to 10. 
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DISCUSSION 


Longitudinal Motion 

Pure longitudinal motion, with no aileron deflection, is shown in 
figure 5 for Doth the conventional and variable -incidence configuration. 
During the pull-up the steady-state value of pitching velocity is the 
same for the two configurations, as it should be, because a lift coef- 
ficient increment of 0.2 was called for d\uring the maneuver. The change 
in steady-state angle of attack is about 7° for the conventional con- 
figuration and very nearly 0*^ when the variable-incidence wing is used. 
The amplitude of the oscillatory motion restating from abrupt deflection 
of the control is much smaller for the variable -incidence -wing 
configuration . 


Combined Motions 


Some preliminary computer runs indicated that for the variable - 
incidence configuration the steady rolling velocity was appreciably 
higher than for the conventional configuration with the same aileron 
deflection. Since a comparison of the two configurations on the basis 
of equal rolling velocities rather than equal aileron deflections is 
more meaningful for a study of the motions, the aileron deflections for 
the variable -incidence configuration were reduced for almost all of the 
runs so as to produce rolling velocities more nearly equal to those for 
the conventional config\nration. 


Time histories.- For the combined longitudinal and lateral motions 
the majority of the calciaations were performed with Cjp = -0.10, which 


corresponds to the airplane with no wing dihedral. These resxats are 
presented in figures 4 and 5 . The aileron deflections were chosen so 
that the largest values used would give roll velocities near the uncou- 
pled yaw frequency (20 radians/sec), according to an estimate from an 
equation for a single-degree-of -freedom roll. For the conventional con- 
figuration with = -0.10 the rolling velocities are only about one- 


half of those from the single -degree -of -freedom calculations. This 
result is primarily due to positive sideslip angles, developed d\nring 
the maneuvers (fig. 4), which produce opposite rolling moments because 
of the dihedral effect. These positive sideslip angles are caused by 
coupling terms other than aerodynamic since On and On were assumed 

Oa T 

to be zero for all calculations herein. (See table II.) Coupling 
effects on the longitudinal motion of the conventional configuration are 
evident in figures 4(b), (c), and (d) by increases in magnitude of q 
and a diiring the time the aileron is deflected. Vfhen the aileron 
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deflection was applied near the first peak in angle of attack (fig. 5) 
the lateral motions of the conventional configuration were more violent, 
with the rolling velocity reversing sign for a short time in one case 
(fig. 5(a)). 

The first calculations with the variable -incidence wing were made 
with the fuselage steady-state angle of attack equal to zero. These 
calculations indicated oscillatory instability for the smallest aileron 
deflection (dotted line, fig. ^(a)) and divergence for the larger aileron 
deflections (dotted lines, figs. 4(b) and 5(a)). The aileron deflections 
for these runs were the same as for the conventional configuration, and 
the higher rolling velocity is evident in figure 4(a) . It is not believed 
that the larger rolling velocities are the cause of the instability, how- 
ever, because as shown later higher rolling velocities were obtained on 
other r\ms without producing instability. Since for the condition of 
zero fuselage angle of attack the airplane principal axis is inclined 
nose down with respect to the flight path, it appeared that poor damping 
of the lateral motion may have been the basic cause of the unsatisfactory 
motion and that only a small amount of coupling produced instability. 

Note in figure 4(a) that the oscillatory lateral motion, which is excited 
by the aileron deflection and then increases in amplitude, subsequently 
damps when the wing deflection is returned to zero. 

In order to obtain better damping of the lateral oscillations the 
fuselage can be trimmed to some positive angle of attack rather than 
zero, as in the first variable -incidence calculations. A constant deflec- 
tion of the horizontal tail will accomplish this trimming; accordingly, 
two additional sets of calculations were made with the fuselage angle of 
attack at I /3 and 2/5^ respectively, of the value of the angle of attack 
of the conventional configuration during the pull-up at a lift coeffi- 
cient of 0.2. Both sets of calculations showed dair^ed motions with no 
divergences, and there was little difference between the two sets; there- 
fore, all motions subsequently discussed for the variable -incidence -wing 
configuration are for the condition of the fuselage steady-state angle 
of attack equal to I /5 of that for the conventional config-uration during 
the pull-up. 

The importance to the coupled motion of the ap term in the side- 
force equation can be seen by observing the effect of the steady-state 
angle of attack on the initial sideslip transients in figures 4(a), 

4(b), and 5(a*) . A decrease in angle of attack from 0.12 radian ( 6 . 9 ^) 
for the conventions.! configuration to 0.04 radian ( 2 . 50 ) for the vari- 
able incidence configuration resulted in large reductions in the initial 
sideslip angles . A further decrease in angle of attack to zero caused 
a reversal in sign of the initial sideslip motions. The angle of attack 
of the principal axis is negative for the latter case. In figures 4(b) 
and 4(a) it is evident that for the zero angle-of-attack case the posi- 
tive rolling velocity and negative sideslip angle are reinforcing each 
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other to the extent that apparent divergences in both the lateral and 
longitudinal motions existed. 

Although the aileron deflections for the variable-incidence con- 
figioration with a positive trim angle of attack are about 25 percent 
smaller than for the conventional configuration^ the rolling velocities 
are slightly higher (figs, k and 5 ) because of the smaller sideslip 
angles developed by the variable -incidence configuration during rolling. 
The yawing velocities during the roll are considerably smaller for the 
variable -incidence airplane. 

For the cases with aileron applied near the first peak in angle of 
attack (fig. 5) the variations in rolling velocity and sideslip angle 
were much less for the variable -incidence configuration (long dash line, 
fig. 5) than for the conventional configuration (solid line, fig. 5) • 

No reversals in direction occurred for the variable -incidence configura- 
tion prior to return of the aileron to neutral. Effects of the lateral 
motion on the longitudinal motion are evident in figures k and 5 ‘the 
variable -incidence configuration, but in no case do these effects appear 
to be larger than for the conventional configuration. 

For the motions in figures 4(a), 4(c), and 5(^) "the amplitudes of 
the lateral oscillations remaining after the ailerons are returned to 
neutral are larger for the variable -incidence configuration than for 
the conventional configuration, while for all other motions illustrated 
the opposite is true. The effect seemed to occur randomly in about l/4 
to 1/5 of the comparable analog runs made for the two configurations 
and is probably a function of the initial conditions existing when the 
ailerons are returned to neutral. 


The large losses in rolling velocity caused by the sideslip angles 
indicated that in terms of airplane flying qualtities the effective dihe- 
dral of the airplane was excessive for the calculations previously dis- 
cussed. Additional calculations were made with reduced to -0.02, 

as might be obtained with negative wing dihedral, and some of these 
results are presented in figures 6 and 7- Note that the p scale in 
figures 6 and 7 has been reduced by a factor of five from that used in 
previous figures. 


The conclusions obtained from the study of motions with 



- 0.10 


are even more strongly evident from the motions with = -0.02. For 

a given aileron deflection smaller losses in rolling velocity caused by 
sideslip angles occurred when was reduced. The variable -incidence 


airplane again produced higher rolling velocities than the conventional 
configuration for the same aileron deflection. With the dihedral reduced 
the variable -incidence airplane had an even greater advantage over the 
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conventional airi^lane and produced maximum sideslip angles only about 
1/5 as large. Tlie yawing velocities were again much smaller for the 
A/ariable-incidence airplane than for the conventional. 

SiiTnmfl-ry plots . - Some of the results of this investigation are sum- 
marized in figures 8 to 10. Figure 8 shows the approximate steady-state 
rolling velocitiess obtained compared with those that would be obtained 
from a s ingle -deg;ree -of -freedom roll. The large losses in rolling veloc- 
ity due to the sideslip developed are evident. For the same aileron 
deflections the rolling velocities for the variable-incidence configinra- 
tion were higher than for the conventional configuration by about 50 per- 
cent for = -0.10 and by about 30 percent for = -0.02. 

The values of maximum sideslip angle shown in figure 9 are signifi- 
cantly less for the variable -incidence configuration than for the con- 
ventional configuration at the same rolling velocity. Although the 
resonance condition was not reached it was closely approached and fig- 
ure 9(b) shows that for the conventional configuration the sideslip angle 
exciarsions were increasing at a rapid rate at the higher rolling veloc- 
ities. The rate of increase with rolling velocity for the variable- 
incidence configuration was much smaller. 

The large reductions in maximum yawing velocity produced by the 
variable-incidence airplane compared with the conventional are illus- 
trated in figure 10. The maximum yawing velocities were not as greatly 
affected by reducing the dihedral as were the maximum sideslip angles. 

Inertia coupling terms . - The terms in the equations listed previ- 
ously which cause coupling of the lateral and longitudinal motions are 
terms proportional to pp, ap, qr, pq, pr, r^, and p2. For the 
airplane configuration and flight conditions used in the calculations 
herein, the results indicated that for equal rolling velocities the 
quantities p, a, r, and the transient portion of q were smaller 
for the variable-incidence-wing type of pitch control than for conven- 
tional tail control. While p2 and the steady-state values of pq 
were the same for the two modes of operation, all other coupling terms 
were smaller for the variable-incidence configuration. This fact was 
confirmed by recorded values of some of these coupling terms obtained 
but not illustrated herein. For a rolling maneuver during a steady 
pull-up (q = constant) the term pq will be the same regardless of the 
type of control used or of the airplane attitude (except insofar as p 
is affected by the choice of axes) . By referring to equation (5) it can 

Iy - ly 

be seen that if q is constant the term 1 pq produces an effec- 

tive negative value of Cn for the type of airplane under consideration 
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larger than as pointed out in reference 3. Negative values 

of Cnp produce destabilizing effects on the lateral motion (refs. 3 

and lo) . Because significant reduction in p and r was obtained 
when the variable -incidence wing was used, the conclusion is reached 
that the destabilizing effect of the pq. term was small when con^ared 
with the beneficial reductions in other coupling terms such as ap . 


Application to Full-Scale Airplane 

As mentioned previously, the flight conditions used in this inves- 
tigation were those existing at a Mach n\imber of 2.0 for the rocket- 
propelled airplane model of reference 8. The motions obtained in this 
study, when nondimensionalized to the forms pb/2V, qc/2V, and rb/2V, 
can be considered as applying to a fxill-scale airplane having dynamic 
similarity; that is, the same value of the relative density factor p 
and the same nondimens ional radii of gyration such as ky^b and ky/S. 

For example, if the model used is assumed to be a l/6-scale model of a 
25,000-pound airplane, then the equivalent airplane altitude is about 
k0,000 feet. The results in figures 3 to 10 may be considered to apply 
directly to such an airplane if the values of the scales for angular 
velocities are divided by 6.56 and the values of the time scale are mul- 
tiplied by the same factor since this factor accoiants for the airplane- 
to-model scale ratio and the variation of speed of sound with altitude. 

For the present investigation the value of Cm^ for the variable- 

incidence -wing configuration was kept the same as for the conventional 
configuration so that the natural frequencies in pitch, which appear in 
frequency ratios that are important in roll coupling (ref. l) , would be 
the same. For a conventional airplane the static stability Chiq^ pro- 
vides stability of the short -period motion and acts as a spring constant 
which must be overcome by the longitudinal control Cmg in order to 

maneuver the airplane. For a variable -incidence airplane the latter 
function disappears, and the criterion for establishing desirable values 
of static stability may be quite different. For example, greater free- 
dom might be possible for selecting the pitch frequency in relation to 
the yaw frequency and thus further reduce the coupling tendencies. 

The condition 0^^ = 0 may be satisfied for any given set of flight 

conditions by the proper choice of size and location of wing and tail, 
as indicated previously . Xt is not to be ejq)ected that ^mg will remain 

zero as the Mach mmiber and altitude change. The value of Cmg might be 

kept relatively small, however, by suitable choice of tail plan form and 
aspect ratio. The effect of such nonzero values of Cmg on the coupled 
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motions wo\ild, of course, be an important factor for investigation when 
assessing the value of variable-incidence-wing control. 


CONCLUDING REMARKS 


The brief investigation described herein indicates that the use of 
a variable -incidence wing as the primary longitudinal control for air- 
planes may offer some advantages over a conventional tail control in 
the matter of reducing undesired dynamic effects during pullout rolling 
maneuvers. Some of these advantages are: smaller pitch transients 

during a symmetrical pull-up; greater rolling velocity for a given 
aileron deflection because of smaller adverse sideslip angles; and a 
reduction in inertia coupling between longitudinal and lateral motions. 
It was indicated that one of the primary effects leading to the reduced 
coupling was the small fuselage angles of attack permitted by the 
variable-wing-incidence type of control. The present study is limited 
in that only one set of flight conditions was used in the study, and 
no attempt to derive the optimum variable -incidence configuration was 
made. It is believed, however, that this investigation indicates suffi- 
cient promise of benefit to be derived from the use of a variable - 
incidence wing that further investigation is warranted. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 5> 195T- 
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TABLE I 

FLIGHT CONDITIONS AND MASS CHARACTERISTICS 

M 2.0 

V, ft/sec 2,200 

q', Ib/sq ft 5,?60 

ft 1,547 

W, lb 500 

l x, slug-ft2 

ly, slug-ft^ 102 

slug-ft^ 104 

1^2^, sliig-ft^ 4 
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TABLE II 


AERODYNAMIC CHARACTERISTICS 


""a 


^^5 


"mo 




- m . 

a 


CZf 

Cz , 


'nr 


' Ht 


"He 


- c 


n: 




Conventional 


Variable incidence 


2.0 

2.0 

0.12 

1.5 

- o . 6 o 

- o . 6 o 

- 0.25 

0 

- 1.10 

- 1.10 

- 0.55 

- 0.55 

- 0.75 

- 0.75 

- 0.10 

- 0.10 

- 0.02 

- 0.02 

0.01 

0.01 

- 0.15 

- 0.15 

0.1 

0.1 

0.2 

0.2 

- 2.0 

- 2.0 

0 

0 

0 

0 
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Figure 1.- Body-axes system. All quantities illustrated in positive 

direction except p . 



Conventional 
Variable incidence 




Figure 2.- Airplane configuration used in study. 
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0 .2 4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 


t , sec 


Figure 5»- Longitudinal response to pitch control deflection. 
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(a) 6a^c = 0.200; &a,v = 0 . 15 ^ and 0.200. 

Figure 4.- Combined lateral and longitudinal motions. Aileron deflection applied at 1.0 sec. 

C, = -0.10. 
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(t)) 5a,c “ O. 5 OO; 6 a, V = 0*250 and O. 5 OO. 
Figure 4.- Continued. 
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5 a,c = 0-^00; = 0.510. 

FigTjre 4.- Continued. 
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rad 
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rad 

sec 
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O. 5 OO; = 0 . 590 . 


Figure 4.- Concluded 
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(a) = 0.500; = 0.231 and O.5OO. 


Figure 5.- Combined lateral and longitudinal motions. Aileron deflection applied after qne- 

fourth cycle of pitch oscillation. C7 = -0.10. 
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Figure 5*- Concluded. 
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Figiire 6.- 


Combined lateral and longitudinal motions. 
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Aileron deflection applied at 1.0 sec. 
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Figure 6.- Concluded. 
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Combined lateral and longitudinal motions. Aileron deflection applied after one-fourth 
cycle of pitch oscillation. Cj^ = -0.02; 6g,^c = 0.400; = O.5o7- 
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Figure 8.- Approximate steady-state rolling velocities 
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(b) Cip = -0.02. 


5.- Mauximum sideslip angles obtained diiring rolling maneuvers. 
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(b) C, = -0.02. 


Figure 10.- Maximum yawing velocities obtained during rolling maneuvers. 
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